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■ Abstract. We present a study of charged lepton mass matrix diagonalization in R- 

, parity violating SUSY. The case in which the bilinear couplings /Lt« have large values is 

' given special attention. 

p; INTRODUCTION 

<d : 

R-parity violating SUSY is a subject which has enjoyed a lot of interest in the past 
few years. The possibility of mixing between particles and superpartners, allowed 
^ ' in these models, makes for very interesting phenomenology. However, due to the 

large number of parameters present, this subject is also quite difficult to study. 
Recently, it has been found (see [1]) that by working in a specific basis (single VEV 
parametrization) the analysis of the fermionic sector is greatly simplified, without 
loss of generality. In this basis, the only non-MSSM parameters that play a role in 
the leptonic phenomenology at tree level are the three RPV bilinear couplings //j. 

The phenomenological consequences of this model in the fermionic sector have 
been extensively studied [2]. This paper aims to detail the technical aspects con- 
cerning the chargino-lepton mass matrix diagonalization. It has been found that for 
large values of the couplings (of order of hundreds GeV and above) this problem 
is not trivial. We should point out that our interest in this issue here is mostly 
theoretical: most of the range of values relevant to the discussion below is not 
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allowed by the experimental constraints (except in special cases; for more details 
on this subject see [2]). 

The framework is the same as in [1]. The mass matrix can be written: 
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where g v = g2V u /V^ = V2M w sinf3, g' v = g2Vd/V% = V2M w cosf3. Here, by m 8 we 
denote the Yukawa masses of the three leptons e, //, r. The physical masses will be 
denoted by m t . 

The aim is to go from the weak interaction fields, in terms of which the lagrangian 
is written initially, to the mass eigenstate (physical) fields, which can be observed 
experimentally. To this purpose, we rotate the left fields by a matrix U L and the 
right fields by a matrix U R . These rotations will also diagonalize the squared mass 
matrices: 

v iA M l jjL = jjR] m r jjR = diag | M 2^ M% 2 ,m%} 

where M L = M\M C , M r = M C M\. 

The mixing between leptons and charginos naturally leads to changes in the 
couplings of these particles to the gauge bosons (see [3]). In the case of Z coupling 
we will have : 
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~A% - 2sin 2 ^, ~A% = 2U* U« + U* U* 



The diagonal (5^) terms are the SM values, while the A terms are consequences 
of the mixing. Being nondiagonal, they lead to anomalous Z couplings and non- 
standard decays {e.g. Z — > e/z, \x — > eee). It turns out that, in most cases, the 
anomalous left coupling is the important one; the anomalous right coupling is pro- 
portional with the product m^m., and very small numerically. In what follows we 
will concentrate on the left rotation matrix. 



I ANALYTIC DIAGONALIZATION OF MASS 

MATRIX 

In this section we will analyze approximate analytical solutions to our diagonal- 
ization problem. The perturbative solution for small //'s has been given in [1]; we 
present it here for comparision with further results: 

Un = Vi t , i = e,n,r (1) 



A = [i M 2 — g v g' v . This formula reveals to us one important fact: the strength 
of the mixing decreases at large tan/3. The phenomenological consequences of this 
behaviour have been analyzed in [1], [2]. 

The fact that the above formula doesn't work at large fi is apparent; indeed, 
as we increase //j, the components Uf{, as given by the above formulae, increase 
indefinitely, which is not allowed by U matrix unitarity. 

It is possible to diagonalize the matrix without assuming that the /ij are small. 
Instead, we will take the Yukawa masses to be small; actually, we will take them 
to be zero in the matrix M L . Then, we get the following solution: 
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where A 2 = A 2 + ct 2/ u 2 , A 2 = A 2 + a 2 (/i 2 + /i 2 ) , A 2 = A 2 + a 2 {[i\ + f4 + A*l) , and 
a 2 = M| + g£ 2 . The first two columns correspond to the two charginos, while the 
last three columns correspond to the leptons (e, //,r respectively). 

Note that, in the approximation used, the physical masses of the particles (eigen- 
values of the matrix M L ) are also zero; as a consequence, the three lepton eigen- 
vectors are degenerate. To get the correct combination, we employ a limiting pro- 
cedure: start from the exact eigenvectors, and let the electron mass, muon mass, 
and tau mass go to zero, in this order. If the lepton mass hierarchy holds also for 
the Yukawa masses (m e << m M << fh T ) it can be shown that this way we get the 
correct result (2). 

Besides lepton eigenvectors, other quantities of interest obtained through diag- 
onalization of the mass matrix are chargino masses. In the small Yukawa mass 
approximation, we get 
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with a 2 = fil+g^+nl (/x 2 , = +/-*!)■ Interpretation of the quantities a\ and a 2 

is straightforward; at large /X5, the mass of the heavier chargino is M Xl ^Jol{ ~ /is, 
while the mass of the lighter chargino is M X2 ps y/a^. Actually, it can be shown 
that the lighter chargino mass increases monotonically from the MSSM value (for 
A*5 — 0) to ^fa^ (for /x 5 — > oo). This behaviour has important phenomenological 
consequences. Consider the fact that the lower limit on the lighter chargino mass 
M\2 > 90 GeV eliminates part of the (M 2 ,/i ) plane in the MSSM. With R-parity 
violating terms, you can expect that this excluded region will shrink; indeed, if we 
make /X5 big enough, it might potentially go away completely. The fact that M\ 2 
is limited above by Ja^ means that some region in the (M 2 , /xq) plane does in fact 
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FIGURE 1. Strength of left-handed Ze/j, anomalous coupling : dotted line: small /i approxima- 
tion (1); solid line: small Yukawa mass approximation (2); dashed line: exact numerical result; 
/ii ratio 1:1:1 



remain excluded, no matter how strong the R-parity is violated. This region is 
given by the equation: 

< 90 GeV 

or, at large tan/5, M 2 < 90 GeV. This result is supported by the exact numerical 
analysis presented in [1], [2]. 

Let's turn back to the lepton eigenvectors and consider the anomalous lepton- 
lepton-Z couplings. For simplicity, let's take the left-handed Ze/i coupling: 

A eil = /ii/i 2 2 (3) 

At small /x's, the strength of the coupling increases with 1^1/^2, while at large //'s it 
decreases like 1 / \Zm§ (see Fig. 1). In between it will reach a maximum value: 

a' 2 

aL _ si jv 

e M max ~ M | + g ,2 

C being a constant which depends only on the ratio /ii//i2- Experimental con- 
straints on anomalous branching ratios or lepton number violating decays (which 



can generally be written in the form A% < Aq) will then be satisfied not only in 
the region of small /i^, but also for large values of 115. Moreover, if in some region 
of parameter space A^ max < A for some process, then that particular process 
won't contribute at all to constraints on /i 5 values. 

The analytic solution (2) derived in this section is not valid at arbitrarily large 
yu's. The approach used to derive the lepton eigenvectors works only as long as 
rh e « fhfj, « m T . Once we get close to the boundaries of the region where the 
diagonalization problem has solutions (see next section) this relation does not hold 
anymore. However, the numerical results for the anomalous couplings Afj show an 
even steeper decrease in this region than that given by (3). As a consequence, for 
sufficiently large /i 5 , constraints from electroweak processes like lepton or Z decay 
dissapear. This region of large 7/5 is excluded on the basis of strong interaction 
processes (ir decay, or neutrinoless beta decay). However, the corresponding case 
involving mainly large /i 3 may be phenomeno logically viable [2]. 



II YUKAWA MASSES 



To be able to perform an exact diagonalization of the mass matrix, we first 
have to find the Yukawa masses of the leptons. Note that this is not the standard 
problem, in which we have a matrix and we have to find eigenvalues. In this case, 
we know three of the eigenvalues (the physical masses m; of the leptons) and we 
have to find some elements of the matrix itself (the Yukawa masses fhi). 

The direct approach would be to use the eigenvalue equations: 

det(M - AJ) = , \ = m 2 e , ml, m 2 T (4) 

This is a three by three system of nonlinear equations, and is not easy to solve even 
numerically (except in particular cases). So, we will use another approach. 

This approach is based on the observation that if we know the solution for some 
values [if , we can find the solution at an neighbouring point fij + dfif. For sim- 
plicity, let's consider the ratio of the /z's fixed, and their magnitude given by a 
parameter t: 

Hi = nVt 

Write the eigenvalue equations: 

E(t,ml,X) = det(M-XI) = 



To an infinitesimal modification in t will correspond an infinitesimal modificaton 
in the Yukawa parameters mf. 

dE dE dmj 

~dt dml ~~df ~ (5j 



Now, we can numerically integrate this system of linear equations from zero to 
whatever value of t we need. 

What about the existence of solutions for our problem? Let's suppose we can 
solve the system (4); in order that the result make sense, we require the solutions 
to be real. This will restrict the allowable values of /ij to some domain D in /ij 
space. What this means, in terms of solutions derived with the help of (5), is that 
we can increase t only as long as we stay inside this domain. When we reach its 
boundary, usually what happens is that the determinant of the system (5) becomes 
zero, and we cannot solve for drhi. 

Another relevant question is if this domain is simply-connected; that is, starting 
from fii = 0, can we reach any point of it with a path formed by connected straight 
lines? In other words, does integrating the system (5) allow us access to all the 
solutions to (4)? 

We do not know the answer to this question for the general case of three leptons. 
But, if we consider the simpler case of only two leptons (presume one of the /ij is 
zero), we can write (4) in the form: 

aifh\ + a2fh\ = s 
a\a2fn\fn\ = p 

with solution: 

(am 2 ) lj2 = S - T \/s 2 -4p 

The quantity A = s 2 — 4p becomes for some value /i 5 = /i 5 max (which gives us 
the boundary for the domain D), and it can be shown that only for /i 5 < /i 5 max 
is A positive (necessary condition for real solutions). Note that /is max is generally 
around a few TeV. 

Another interesting issue is the problem of lepton mass hierarchy in this model. 
In the Standard Model (or MSSM) we have m e << << fh T . These relations 
need not hold in our R-parity violating scenario. Take, for example, the two lepton 
mixing case presented above. If \i\ = /i 2 , then, at /is = /is max , we have fh\ = m 2 
(this will happen at quite large /is values, though; /i 5 max ~ 7712/1711). If \i\ > /i 2 , it 
is even possible to get rh\ greater than m 2 . The next question is if this behaviour 
holds for the general case of three lepton mixing. The possibility of finding points in 
parameter space where the three Yukawa masses are of the same order of magnitude 
(or maybe even equal) is particularly interesting. Unfortunately, the technical 
difficulties encountered in working with the nonlinear system (4) have stopped us 
from getting an answer to this question so far. 

Ill CONCLUSIONS 

For large R-parity violating terms, the mixing between charginos and charged lep- 
tons has different characteristics than at small /ij. We have presented approximate 



analytical expressions for both regimes, which can help to understand numerical 
results. 

The problem of existence of solutions of the system (4) - finding the Yukawa 
masses so that three of the mass matrix eigenvalues will be equal to the physical 
lepton masses - is still unsolved for the general case of three lepton mixing (although 
within the phenomenologically viable region, numerical solutions are always suc- 
cessfully obtained in [2]). For two lepton mixing it can be solved, and it has been 
shown that there are regions in parameter space where the Yukawa masses of the 
two particles are of the same order of magnitude. 
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